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Abstract 

The interaction between vortex density waves and high-frequency second sound in coun- 
terflow superfluid turbulence is examined, incorporating diffusive and elastic contributions 
of the vortex tangle. The analysis is based on a set of evolution equations for the en- 
ergy density, the heat flux, the vortex line density, and the vortex flux, the latter being 
considered here as an independent variable, in contrast to previous works. The latter fea- 
ture is crucial in the transition from diffusive to propagative behavior of vortex density 
perturbations, which is necessary to interpret the details of high-frequency second sound. 

1 Introduction 

The study of vortex tangles in superfluid helium II has received much attention during the 
last decades, both because of its own interest and as a first step to understand the classical 
turbulence. Many theoretical and experimental studies on superfluid arguments have enhanced 
the knowledge on these intricate phenomena: experimental studies have allowed direct results 
and confirmed theoretical intuition, while, on the other hand, theoretical studies are important 
not only as a guide for experiments but also as an explanation of the experimental results [T]- 

It is known that the presence of a heat flow in superfluid helium II causes the formation 
of quantized vortex lines, which move inside the superfluid until a stationary situation is 
reached, and whose presence is usually investigated by second sound waves pQ-[3]. But, these 
waves interact with the vortex lines, causing not only a modification of their velocity and an 
attenuation of it, but also a modification of the vortex lines profile and of their motion. In 
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the last years, the study of non-stationary and inhomogeneous turbulent states has attracted 
much attention [4]-[7j. The vortex lines and their evolution are investigated by second sound 
waves, so that it is necessary to analyze in depth their mutual interactions. In particular, 
high-frequency second sound may be of special interest to probe small length scales in the 
tangle, which is necessary in order to explore, for instance, the statistical properties of the 
vortex loops of several sizes. In fact, the reduction of the size of space averaging is one of the 
active frontiers in second sound techniques applied to turbulence, but at high-frequencies, the 
response of the tangle to the second sound is expected to be qualitatively different than at low 
frequencies, as its perturbations may change from diffusive to propagative behavior [2|-[9]. 

In a previous paper [3] a thermodynamical model of inhomogeneous superfiuid turbulence 
was built up with the aim to study the mutual interactions between second sound and the 
vortex tangle. The fundamental fields of the model were the density p, the velocity v, the 
internal energy density E, the heat flux q and the average vortex line length per unit volume 
L. In a successive paper [9], starting from this model, a semiquantitative expression for the 
vortex diffusion coefficient was obtained and the interaction between second sound and the 
tangle in the high-frequency regime was studied. In both these works, for the sake of simplicity, 
the diffusion flux of vortices J was considered as a dependent variable, collinear with the heat 
flux q, which is proportional to the counterflow velocity V, defined as V = v n — v s , v n and 
v s being the velocities of the normal and superfiuid components, respectively. 

But, in general, this feature is not strictly verified because the vortices move with a velocity 
V£, which is not collinear with the counterflow velocity. Indeed, using the vortex filaments 
model proposed by Schwarz in |10j-|12|. which describes the vortex line by a vectorial function 
s (£>^)> £ being the arc-length measured along the curve of the vortex filament, the velocity of 
the vortex element is pQ-[3] 

v L = v si + as' x (V - Vi ) - q's' x [s' x (V - v ; )], (1.1) 

where a and a' are temperature-dependent friction coefficients between the normal fluid and 
the vortex line, s' the unit vector tangent along the vortex line at a given point, and v s ; = v s +Vi 
the "local superfiuid velocity", sum of the superfiuid velocity at large distance from any vortex 
line and of the " self- induced velocity" , a flow due to all the other vortices including other parts 
of the same vortex, induced by the curvature of all these lines. In (11. 1|) the prime indicates 
the derivative with respect to the arc- length £, that is s' = ds/d^. In the "local induction 
approximation", the self-induced velocity Vi is approximated by pQ-[3] 

v^/^xs'W with /3 = ^ln(^), (1.2) 

with c a constant of the order of unity, k the quantum of vorticity, given by k = h/m, h Planck's 
constant and m the mass of a helium atom, ao the size of the vortex core and s" = d 2 s/d^ 2 
the curvature vector. The intensity of Vi is |vi| = fi/R, with R the curvature radius of the 
vortex line. The coefficient /3 is linked to the internal energy per unit length of the vortex line 
(the tension of the vortex line) by the relation ty = p s n/3. 

Using the local induction approximation, equation (jl.ip can be written as 

v L = (1 - a')(3s' x s" + a/3s" + v s + as'xV- a's' x (s' x V). (1.3) 

Consider now a mesoscopic portion of turbulent superfiuid, contained in a small volume A, 
which contains a small vortex tangle. In the following, the vortex velocity < vl >, averaged 
in A, will be denoted with v tangle = v tan . 
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Integrating over the volume A, recalling that in counterflow situations, i.e. for p n v n + 
p s v s = 0, it is v s = — (p n /p) V, and supposing that in A the counterflow velocity is constant, 
one gets 

vtan = (1 - a')j3 < a' x s" > +a(3 < s" > -— V + a < s' > XV + a' < U - s's' > -V, (1.4) 

P 

where < • > stands for the averaged value of the vector in A. 

Suppose that in the small volume A the vortex tangle is homogeneous. If the considered 
volume A is sufficiently far from the walls, and therefore does not contain pinned vortices, one 
can suppose < s' >= and < s" >= 0, obtaining 

v ta „ =< v L >= -^V + la'n s • V + (1 - a')/3Ic u (1.5) 
P 3 

where the tensor II s = § < U — s's' > was studied in [13], while the vector I = ^ig/fj^ and the 

scalar c\ = ^ L3 / 2 were introduced by Schwarz [12j in a microscopic approach to the dynamics 
of superfluid vortex tangles. From f)1.5|) one sees that vt an || V only when II s = U, i.e. when 
the distribution of s' in the tangle is isotropic, and the anisotropy vector I is collinear to 
the counterflow velocity V. However, experiments and numerical simulations show that these 
hypotheses are only approximately verified. 

In the hypothesis of isotropy of the tangle (II s = U), the assumption v tcm = implies 

1 || V. But, in general, one could have vt an 7^ and I not collinear to V, which means 
vtan = AV + B1, with A and B suitable coefficients coming from the relation (jl.5j) . Therefore, 
the hypothesis vt an collinear with V is not in exact agreement with the microscopic results of 
the vortex filament model. For this reason, the aim of this paper is to build up a model of 
inhomogeneous counterflow superfluid turbulence, in which the flux of the vortex lines, which 
is parallel to vt an , is taken as an additional independent variable. 

In Section 2 a relaxational generalization of the diffusion equation for vortex lines is pro- 
posed, analogous to the well known Maxwell-Cattaneo generalization of Fick, Fourier, Ohm 
and Newton-Stokes classical transport laws |144 1 15 j . There, the corresponding generalization 
of the entropy in order to achieve compatibility between the relaxational transport laws and 
the second law of thermodynamics in a simplified model is studied, in which only the line 
density L and its diffusion flux J are fundamental variables. This simplified model allows us 
to understand the contribution of J to the nonequilibrium entropy but does not describe the 
interaction between vortex density waves and second sound. For this reason, in Section 3 a 
more general analysis, including as independent variables of the theory not only the line den- 
sity L and the diffusion flux of vortices, but also the energy and the heat flux is undertaken. 
The mathematical formalism is physically motivated to explore the interactions between heat 
waves and vortex-density waves. In Section 4, the physical meaning of the several coefficients 
appearing in the model are examined, and in Section 5 wave propagation in this more complete 
model in uncoupled and coupled situations is studied and the results are compared with those 
obtained in [HE]. 

2 Simple formulation of vortex- density waves 

This Section aims to provide a simplified view of the behavior of turbulent vortices in inho- 
mogeneous situations, with special emphasis on the transition from diffusive behavior at low 
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frequencies to propagative behavior at high frequencies. This example provides the physical 
motivation for the wider treatment presented in Section 3. 

In this Section the constitutive equation for the diffusive flux of vortex lines is generalized, 
by including relaxational effects due to their inertia. Thus, one will be interested in the 
evolution of inhomogeneous vortex tangles or in the response to external perturbations inducing 
local changes in the vortex line density. Usually, in counterflow situation (i.e. under vanishing 
barycentric velocity) one considers homogeneous vortex tangles and the evolution equation of 
L is assumed to be the well-known Vinen's equation [I]-[3], [16] 

^ = a v VL 3 / 2 - ? v kL 2 = a L , (2.1) 

where ai stands for the net vortex production per unit volume and time and a v and f3 v are 
numerical coefficients. 

When one assumes inhomogeneous vortex tangles, with L differing from point to point, a 
further contribution should be added to (|2.ip . thus becoming 

§— VJ + «. (2.2) 

In (|2.2p . it appears a new quantity, J, the flux of vortex lines. In principle, one could expect 
that J will be related to the gradient of L, in analogy with the well-known Fick's law for matter 
diffusion. In fact, in [3] a thermodynamic formalism leading to such a transport equation is 
studied. In some occasions, however, inertial effects may be relevant, and they will contribute 
to the constitutive law for J with a relaxation term. 

The aim of this Section is to write an evolution equation for J which is compatible with 
the second law of thermodynamics. The dependence of the constitutive relations on the tem- 
perature T and on the heat flux q will be neglected. In the following Section this simplified 
hypothesis will be abandoned. To achieve a consistent thermodynamic formalism, with a pos- 
itive definite entropy production, the entropy must be extended by including J in the set of 
independent variables, as in extended thermodynamics. The corresponding generalized Gibbs 
equation is [H] 

pds = -T~ 1 p L dL - T~ l aJ ■ dJ, (2.3) 

where s is the entropy per unit mass, a a coefficient that will be identified below, and p L the 
vortex chemical potential. Equation (|2.3[) can be written in terms of time derivatives as 

pis = -T~ 1 p L L - T~ l a3 ■ j. (2.4) 
Substituting ([221 in ([23]) one has 



ps + V- (-T~V L J) = J 



V(T~V L ) - T -1 dJ -T~ L p L a L , (2.5) 



i 



where —T~ 1 p L J may be interpreted as the vortex contribution to the entropy flux, and the 
term on the right-hand side as the entropy production. To ensure the positive character of the 
diffusion contribution to the latter, one writes 
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V(T~V L ) - T- X a3 , (2.6) 
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with 7 a positive phenomenological coefficient. A relaxation time tj may be denned as tj = 
7T~ 1 a. With this identification of the relaxation time, (|2.6p may be written as 

Tjj + J = -DVL, (2.7) 

with D a vortex diffusion coefficient identified as D = ^T^ 1 ^^. By using dimensional analysis, 
tj can be taken of the form tj = (71KL) -1 , where 71 is a dimensionless phenomenological 
coefficient. For fast variations of J — namely, in high-frequency experiments — the first term 
on the left-hand side of (j2.7|) becomes dominant. 

Combination of (|2.7p with (|2.2p . neglecting second-order terms in VL, yields 

tjL + L = DV 2 L + Tja L + a L . (2.8) 

For low values of the frequency, the first term of the equation (|2.7p is negligible and one has 
a reaction-diffusion equation, whereas if the frequency is high, the first term is dominant and 
(12. 8p yields a wave equation for L. 

This brief presentation, based on the simplest version of the so-called Extended Thermo- 
dynamics [H], has been enough to give us a simple introduction to the topic studied below, 
namely, the transition from a diffusive behavior to an undulatory behavior of the vortex tangle 
and to the new physical features related with J as an independent variable. However, the evo- 
lution equation for J is introduced as an additional equation to the previous system proposed 
in [1], ignoring the couplings with other possible phenomena, mainly, high-frequency second 
sound. Indeed, in this case, the vortex tangle will not behave as a diffusive system but as an 
elastic matrix, and the dispersion relation for second sound will be changed with respect to 
the previous model studied in [1]. It is need to know in detail these changes in the dispersion 
relation, because of the instrumental importance of second sound as a probe for the vortex 
tangle. 



3 Balance equations and constitutive theory 

To deal with sufficient generality with the interactions between second sound and the dynamics 
of the vortex tangle, one builds up a thermodynamical model of inhomogeneous counterflow 
superfluid turbulence, which chooses as fundamental fields the energy density E, the heat flux 
q, the averaged vortex line length per unit volume L, and the vortex diffusion flux J. Because 
experiments in counterflow superfluid turbulence in the linear regime are characterized by a 
zero value of the barycentric velocity v, in this paper one does not consider v as independent 
variable. In a more complete model v, and p will be also fundamental fields. 

It is known that the heat flux q is linked to the counterflow V through the relation 
q = p s Ts\ , and here one prefers choosing q as variable because it is the macroscopic variable 
appearing in the balance equation for the energy, and it may be controlled in the experiments. 

Consider the following balance equations 

d t E + d k q k = 

d t qi + d k Jf k = a\ (3 ^ 

d t L + d k Jk = <?L 
d t Ji + d k Fik = of 

where dt stands for d/dt and dk for d/dxk, E is the specific energy per unit volume of the 
superfluid component plus the normal component plus the vortex lines, Jf- the flux of the heat 
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flux, Jj the flux of vortex lines, and F{j the flux of the flux of vortex lines; cr|, &L and af 
are the respective production terms. Since in this work one is interested to study the linear 
propagation of the second sound and vortex waves, the convective terms have been neglected. 

If one supposes that the fluid is isotropic, the constitutive equations for the fluxes and 
Fij, to the first order in qi and Jj, can be expressed in the form 

Jf k = P(E,L)S ik , 

F ik = i;(E,L)S ik . { °- ' 

Restrictions on these relations are obtained imposing the validity of the second law of ther- 
modynamics, applying Liu's procedure [151 El- This method requires the existence of a scalar 
function S and a vector function J k of the fundamental fields, namely the entropy per unit 
volume and the entropy flux per unit volume respectively, such that the following inequality 

d t s + d k jjf - a e [d t E + d kqk ] - Af [d tqi + d k j q ik - af] 

(3.3) 

-A L [d t L + d k J k - a L ] - Af [d t Ji + d k F ik - af] > 0, 

is satisfied for arbitrary fields E, qi, L and Jj. In this inequality, which expresses the second 
law of thermodynamics, S and J k are objective functions of the fundamental fields. In order to 
make the theory internally consistent, one must consider for S and J k approximate constitutive 
relations to second order in qi and Jj 

S = S (E,L)+S 1 (E,L)q 2 +S 2 (E,L)J 2 +S 3 (E,L)q l J l , J s k = ^(E, L)q k +<f> J (E, L)J k . (3.4) 

The quantities A E , A|, A L and Aj are Lagrange multipliers, which are also objective 
functions of E, qi, L and Jf, in particular, one puts 

A E = A E {E, L, q u Ji) = A E (E, L) + Af(E, L)q 2 + A E (E, L)J 2 + A E (E, L) qi J u 

A L = A L (E, L, q h Ji) = A L (E, L) + A L (E, L)q 2 + A L (E, L)J 2 + A L (E, L) qi J h 



A] = \\iqi + A12 Ji and Af = \ 2 iqi + A 2 2</i, (3.5) 

with X mn = \ mn (E,L). The constitutive theory is obtained imposing in (I3.3P that the coeffi- 
cients of all derivatives vanish. Imposing that the coefficients of the time derivatives are zero, 
one obtains 

dS = A E dE + A q l dq i + A L dL + A J i dJ i . (3.6) 

Note that S = ps, therefore this equation generalizes equation (|2.3p when energy and heat 
flux variations are taken into account. In the same way, imposing that the coefficients of space 
derivatives vanish, one finds 

dJl = A E dq k + A q dJ? k + A L dJ k + A{dF lk . (3.7) 
Substituting now ([321), (J33J) and ([33]) in (|3.6ti3.7p . one gets 

S\ = -\i\, #2 = 2^22) S3 = A12 = A21, (3.8) 
^ = A E , <P J = A L , (3.9) 



6 



dSo 



AndE + AndL, 



dS^ 



Af dE + AfdL, 



(3.10) 



dS 2 = AfdE + AndL, dS 3 = A^dE + AndL, 



(3.11) 



d(p q = \ n df3 + Midi), dcp J = X 12 df3 + \ 22 dip. (3.12) 

In particular, one obtains to the second order in q and J the following expressions the 
entropy and for the entropy flux 

S = S + ^X u q 2 + \\22 J 2 + Ai29i Ji, J s k = A$q k + A%J k . (3.13) 
It remains the following residual inequality for the entropy production 



a 



s - Afaf + A L a L + Afa{ > 0. (3.14) 



Now, the relations obtained are analyzed in detail. One first introduces a generalized 
temperature as the reciprocal of the first-order part of the Lagrange multiplier of the energy 



dSo 
dE 



(3-15) 



Observe that, in the laminar regime (when L = 0), Aq' reduces to the inverse of the absolute 
temperature of thermostatics. In the presence of a vortex tangle the quantity (|3.15|) depends 
also on the line density L. 

As in QJ. writing equation (|3.10P i as 

dS = ^dE + Aq <±L = AdE - ^dL, (3.16) 

one can identify the quantity — Aq/Aq 5 = —TAq with the chemical potential of vortex lines 
(near equilibrium) 

- TAq = n L . (3.17) 
From (|3.16p one obtains the integrability condition 

dE _ _2 9 / „J 



Neglecting in (|3.6p second order terms in q and J, and using relations (|3.8p . (|3.15p and 
(|3.17p . the following expression for the entropy density S is obtained 

1 fi L 

dS = —dE - —dL + Xnqidqi + A 22 JjdJ, + \\ 2 {Jidqi + qidJi). (3.19) 
Consider now equations (|3.12p . which one rewrites using (|3.9p and (|3. 1T[) as 

d (l) = X n dP + A 2 i#, d (—) = XndP + \ 22 d^. (3.20) 



From these equations, one obtains the following relations 
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df3 



"2 



d = + 



U2 



A11A22 — A12 2 A11A22 — A12 2 V T 



dip 



V12 



rr/ [ 4 



Ai 1 A90 — Aio 2 V T 



from which, putting 



one obtains 



d/3 9/3 _ 



A11A22 — A12 2 \ T J 



dtp dip 

df = r1, ~WL = V ' 



(3.21) 
(3.22) 

(3.23) 



A12C + A22?/ 



A11X + A2i^ = 0, 
d ( fJ-o 



dT 



and also 



Ane + A 2 i?? = -^2, (3.24) 
Ai2X + A 22 ^ = 4r (-4-) > ( 3 - 25 ) 

(3.26) 

X = fWdL' U = ~fW~dL' (3 ' 27) 

where A^ = A11A22 — A12 2 , and v is a positive coefficient because it is the square of the velocity 
of the vortex wave, as it will be shown in the next section. 
Finally, one obtains for the entropy flux 



1 A +A 9 
-7^22 + Ai2^ I — 



1 



1 A12 drf 



7fiM2 ~ An— I J 
1 An di4 



J s - l n ^ T 



(3.28) 



which is analogous to the usual expression of the entropy flux in the presence of a mass flux and 
heat flux, but with the second term related to vortex transport rather than to mass transport. 

Observe that the expression of the entropy flux (|3.28p obtained in this Section is in agree- 
ment with (|2.5p when the dependence on the heat flux is neglected. In the same way, comparing 
the expression of the entropy (|3.19p with the generalized Gibbs equation (|2,3p . proposed in 
the simplified model in Section 2, and keeping in mind equation (|2.6p one gets 



A 



22 



-aT 



TJ 

7 ' 



thus furnishing a physical meaning of the coefficient A22 appearing in previous equations. 

Finally, substituting the constitutive equations (|3.2p in system (|3.ip , and using the relations 
(|3, 21113, 2T|) . the following system of field equations is obtained 
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(3.29) 



' d t E + d jqj = 
dtQi + &jT + X diL 
dtL + djJj = a L 
{ d t J i + r ] d t T + ud i L = af. 

The coefficients 7 and r/ describe cross effects linking the dynamics of q and J with L and T, 
respectively. Thus, they are expected to settle an interaction between heat waves and vortex 
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waves, whose study is one of the aims of the present work. The production terms a must 
also be specified. Regarding of, since only counterflow situation is considering , a simplified 
expression, already noted in literature [I8j LTH] , is assumed 

a* = --kB hv LU s ■ q, (3.30) 
3 

where Bhv = —ol is the Hall-Vinen coefficient [1] and II s = f < U — s's' > is the symmetric 
tensor mentioned in Section 1. If one assumes isotropy in the plane yz, this tensor II s can be 
written as [131 [20] 

/2a \ 
IT = - 1 - a , (3.31) 
2 V l-o / 

where < a < | is a parameter characterizing the anisotropy of the tangle such that < 
s'y >=< s' 2 >= a and < s' x 2 >= 1 — 2a. If the tangle is completely anisotropic, as in the case 
of a regular array produced by rotation, then a = 0, whereas if it is isotropic then a = g. This 
term describes a friction force when q is orthogonal to vortex lines and null force when it is 
parallel to them. A further dissipative term, proportional to the binormal vector I introduced 
in (|1.5p . could be added to right-hand side in (|3.30p (see the Appendix of Ref. [1]), but it is 
neglected here because its contribute is small compared to the right-hand side of (13.301) . For 
the production term a L , one chooses the Vinen's production and destruction terms, equation 
(|2.ip . which one can write in terms of L and of the absolute value of q using the relation 
q = p s TsV 

a L = -BL 2 + AqL 3/2 , (3.32) 

where A = a v /p s Ts and B = /3 v k. For the production term of vortex line diffusion, one 
assumes the following relaxational expression (see the relation below equation (|2.7[) ) 

a J = - 7i kLJ = — , (3.33) 
tj 

where the positive coefficient 71 can depend on the temperature T; with this expression, in 
isothermal situations, one would have a diffusion coefficient given by D = tjv = — ^ . 
Note that in (|3.30p and (|3.33p one has assumed that the respective production terms of q and 
J depend on q and J, respectively, but not on both variables. In more general terms, one 
could assume that both production terms depend on the two fields q and J simultaneously. 

Analyzes, now, the entropy production (I3.14p which, with the expressions of the production 
terms defined above, becomes 

a s = -\ n w (IT • q)< Qi ~ A12 (ti? (IT • q). + J { - \ 22 JikJ? ~ jr (~BL + AqL 1 ' 2 ^ > 0, 

(3.34) 

where w = ^kBhv an d 7 is the positive phenomenological constant defined in (12. 6p . Looking 
at the expression (I3.34p . one notes that the entropy production a s is positive when a suitable 
choice of the coefficients An and A12 is made. In particular, assuming that II s = U, a s in 
(|3.34p is a quadratic form on the variables \q\, L 1 ! 2 and \ J\. Therefore, (|3.34p is verified if the 
matrix 



9 



-t^h l A ±fi L B ' j (3.35) 

-|Ai 2 (k7 + 7i«) -A 2 27i K 

is semidefinite positive. This implies that the coefficient An has to be negative and, being 
d^/dL > 0, from the relation (EHD) one deduces that N is positive. 
Observe that the field equations for q and J can be written also as 



""' A22 ^ + (£) = -1-BhvL [(3a - l) qi 5u + (1 - a) Qi ] , (3.36) 



dt NT 2 dxi N dxi \ T r 

dJ, A 12 dT An d ffx L \ 

^ + NT^ l -^r^vf)= " 71KLJ - (3 - 37) 

Comparing the equation (|3.37p with (|2.6p . one deduces that they can be identified with each 
other if one puts Ai 2 = 0, A 2 i = 0, A 22 = — T~ 1 a and 7]" 1 = T~ l ^anL. Observe also that 
under this hypothesis equation ()3.36|) becomes 



^ + H— = a 

dt dxj 



where ^ = — x ^ T ~i ■ This latter equation is identical to that used in Refs. [TBI H9] , where the 
fields L and J{ were considered as dependent variables. 

4 Physical meaning of the coefficients of proposed equations 

In order to determine the physical meaning of the coefficients appearing in equations (|3.29p ~ 
(|3.33|) . concentrate first the attention on the equations for L and J, 

3 t L + d l J l = a L (4.1) 
d t Ji + r]d t T + udiL = a Jl = -71^. (4.2) 
Supposing that J varies very slowly, Eq. (j4.2j) gets the form 

/; VT V — VL. (4.3) 



71 kL 71 kL 

Substituting it in (|4.ip . one obtains 

d t L = — T1 — V 2 T+ — V —V 2 L + a L . (4.4) 
71 kL 71 kL 

It is then seen that the coefficient ^ = D% represents the diffusion coefficient of vortices 
as already introduced in (|2.8p in a simpler setting. Coefficient — ^ = D 2 may be interpreted 
as a thermodiffusion coefficient of vortices because it links the temperature gradient to vortex 
diffusion. In other terms, this implies a drift of the vortex tangle. Detailed measurements 
have indeed shown [1, pag.216] a slow drift of the tangle towards the heater; this indicates 
that 77 < and small. The hypothesis 77 = corresponds to D 2 = 0, i.e. the vortices do not 
diffuse in response to a temperature gradient. Now, focus the attention on the equations of T 
and q 
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pc v d t T + pe L d t L + d iqi = 0, (4.5) 
d t q l + tdiT + x diL = <r q = -^kB hv LU s • q. (4.6) 



Supposing dtqi negligible in (|4.6p . one gets 



that is 



(IT • q), = ^— VT ^^VL (4.7) 

kBhvL kBhvL 



3£ /TTS\-1V7T ^ i'TTS\-l 



The first term in (14. 8p may be identified as a tensorial thermal diffusivity, and the second 
one is analogous to Soret diffusion term, which describes a coupling between heat flux and 
concentration gradient in usual fluids mixtures; here, instead of the concentration of a chemical 
species, one has a vortex density gradient. 
Substituting (14.80 in ([13]) , one gets 



pc v d t T + pe L d t L = (IIT 1 V 2 T - (IT)" 1 V 2 L, (4.9) 



and assuming isotropy one gets 



d,T= — 



pc v 



7iftL kBhyL 



v 2 t + — 



7iAtL kBhyL 



V 2 L, (4.10) 



From the relations (|4.4p and (|4.10p . and from the positive character of the vortex diffusion 
coefficient and of thermal conductivity one deduces 

— >0 (4.11) 

7i 

and 

^V^<0. (4.12) 
kBhvL "fl^L 

Thus, despite of the high number of coefficients appearing in (|3.29p - (|3,33p . one has been 
able to provide a physical interpretation for many of them, which would allow for their respec- 
tive measurements in suitable experiments. 



5 Interaction of second sound and vortex density waves 

In this Section wave propagation in counterflow vortex tangles is studied, with the aim to 
discuss the physical effects of the interaction between high-frequency second sound and vortex 
waves. Expressing the energy E in terms of T and L, the system (|3.29p becomes 

pc v d t T + pe L d t L + djqj = 0, 
fa + &T + xdiL = *l 

8 t L + d j J J =a L , { ' 

d t J t + rAT + vd t L = a J \ 

where cy = 8tE is the specific heat at constant volume and ej, = d^E. These equations are 
analogous to those proposed in [4] except for the choice of Jf. in fact here Jj is assumed to 
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be an independent field whereas in [4] Jj was assumed as dependent on q. L . However, at high 
frequency, Ji will become dominant and will play a relevant role, as shown in the following. 
A stationary solution of the system (15. ip . with the expressions of the production terms 

(EDBal) , is 



TnU) 



qo = (901,0,0) 

kB hv 



L = Lv, 



A 2 



Qai, 



L aqoiXi, J c 



B 2 
( kB hv 



aq 01 ,0,0 



(5.2) 
(5.3) 



£ V £71 « 

with q i > 0. 

The quantities (13.30f) . (13.32H and (|3.33p can be approximated around the stationary solu- 
tions in the following way 



1 



kB hv [(3a - l)aa + (1 - a)] (g i0 (L - L ) + L qi 



2BLn 



3 AT 1 / 2 

2 AL o 9oi 



(L - L Q ) + ALo /2 q • (q - q 



and 



-7i/cL J - 7ik(L - L )J , 



(5.4) 
(5.5) 

(5.6) 



where the subscript denotes the stationary values for q, L and J. 

Now, consider the propagation of harmonic plane waves of the four fields of the equation 
(|5.ip in the following form 



T = T (x) + fe i( - Kn -*-^ 
q = q + q e i(Kn- x -ut) 
L = L + l e i(K^-wt) 

J = J + 3e^ Kn ^~ wt \ 



(5.7) 



where K = k r + ik s is the wave number, u; the real frequency, n the unit vector along the 
direction of the wave propagation, and the oversigned quantities denote small amplitudes of 
the fields, whose product can be neglected. 

Substituting (|5.7p in the system (|5.ip , the following equations for the small amplitudes are 
obtained 



-uj[pc v ] T 



Ld[pe L ] L + Kq - n = 
\kB hv Lq ((3a - l)cici + 1 - a)] q + ^KTn 
- (-xoKn + iaKBnvCLo) L = 



-qj - 1 [2BLq - ^aLq q 01 ^ 
(— oj — fyiKLo) J + tjqKiiT + {vqKiy — i^inio) L = 

where ci is their unit vector along the first axis x\ and ciCi is the dyadic product. Note that 
the subscript refers to the unperturbed state; in what follows, this subscript will be dropped 
out to simplify the notation. 



\AL 



V 2 , 



L + KJ ■ n + iAL^qx 



(5.8) 
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First case: n parallel to q 

Now, impose the condition that the direction of the wave propagation n is parallel to the heat 
flux qo, namely n = (1,0,0). Through these conditions the system (I5.8P becomes 



—upcyT + Kq\ — bjpeiL - 
^KT — (oj + ianBuvL) q"i 







iAL 3 / 2 qi - (u + iT L x ) L + K J x 
r)KT + [yK - i-fiKJi) L + (-w - 



= 

ryiz-cL) J i 



^kB hv L(1 
«7ikL) J 2 = 
?7ikL) J 3 = 



(5.9) 








where 



-AL 1 ' 2 q 1 
2 y 



r L 1 = (^2i?L 

Note that the transversal modes, those corresponding to the four latter equations, evolve inde- 
pendently with respect to the longitudinal ones, corresponding to the four former equations. 

One will limit the study to the case in which u and the modulus of the wave number K 
assume values high enough to make considerable simplification in the system. Indeed, it is 
for high values of the frequency that the wave behavior of the vortex tangle can be evidenced 
because the first term in (|5.1b ) will become relevant, as shown in Section 2. Note that the 
assumption \K\ = \k r + ik s \ large refers to a large value of its real part k r , which is related to 
the speed of the vortex wave, whereas the imaginary part k s , corresponding to the attenuation 
factor of the wave, will be assumed small. 

This problem is studied into two steps: first assuming \K\ and u extremely high to neglect 
all terms which do not depend on them. Then, the solution so obtained is perturbed in 
order to evaluate the influence of the neglected terms on the velocity and the attenuation of 
high-frequency waves. 

Step I: Under the mentioned assumptions the system (|5.9I) becomes 



—ujpcyT + k r q~\ — ujpe^L 
£k r T — uqi + xk T L = 
— ojL + k r J\ = 
rjk r T + vk r L — ujJ\ = 







(5.10) 



—u>q2 = 0, — ujqs = 0, —L0J2 = 0, — oj J3 = 
Denoting with w = u/k r the speed of the wave, the following dispersion relation is obtained 



w 



^ 1/2 

where V2 = (—XnT pcv) is the second sound speed in the absence of vortex tangle 
[H[T8l[T9] and from (|3.24b ) it is related to the coefficient £ by the relation £ = V^pcy— Xizv/Xn. 
Further, if one assumes that the coefficient rj is zero 

,.0 * * ( , 12) 



V 
pc v 



peL 



w 2 + V%v 



0, 



(5.11) 
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then the dispersion relation (15. llf) has the solutions 



Wx,2 = ±V2, W 3A = ±Vv, 

to which correspond the propagation modes shown in Table [TJ 



(5.13) 



L -- 
Ji 





= 



T ■- 

Qi 
L -- 
Ji 



l 

pc v 



Table 1: Modes corresponding to second sound velocity and vortex waves, respectively. 



As one sees from the first column of Table[H under the hypothesis (|5.12p the high-frequency 
wave of velocity wx,2 = ±V 2 is a temperature wave (i.e. the second sound) in which the two 
quantities L and J\ are zero, whereas in the second column the high-frequency wave of velocity 
11^4 = iby^ is a wave in which all fields vibrate. The latter result is logic because when the 
vortex wave is propagated in the superfluid helium, temperature T and heat flux q\ cannot 
remain constant. This behavior is different from that obtained in [9], because using that model 
in the second sound also the line density L vibrates. In fact, there the flux of vortices J was 
chosen proportional to q, so that vibrations in the heat flux (second sound) produce vibrations 
in the vortex tangle. Experiments on high-frequency second sound are needed to confirm this 
new result. 

Step II: Suppose that the terms of the system (15. 9h . which don't appear in the system 
(|5.10j) . and the term r\ are small enough to be considered as perturbations of the velocity w 
of the wave and of the attenuation term k s of the wave number K. Substituting the following 
assumptions 

w = — = w + 5 and K = k r + ik s 



in the system (|5.9p . one find the expression (|5.13p . at the zeroth order in 5 and k s , whereas at 
the first order in 5 and k s , one obtains 



Itfl 2 



2pc v wU - w% A 



X 



w 



3,4 



(5.14) 



W 3 A 



and 



1 + 



ki 1 ^ 



2pcy [ w l t 2 ~ W 3,' 



PtL 



X 



w 



3,4 



aKLB HV AL z I 2 {w\ 2 pe L - x) 



2w x ,2 

, (3,4) _ ^ L 7i + Tj 
2w 3 ,4 



+ 



2(<2 



w. 



3,4 



^3,4, 



1 AL 3 / 2 (wl 2 pe L - X ) , JxK-n. 



2 ( w* 



w 



3.1 



+ 



2^,4 



(5.15) 

(5.16) 
(5.17) 
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Observe that in this approximation all thermodynamical fields vibrate simultaneously and the 
attenuation coefficients k s are influenced by the choice of J as independent variable, as one 
easily sees by comparing expressions (15,16fJ5~T71) with those obtained in [9\. Looking at these 
results, in particular the two speeds (I5.14I45TT51) . one sees that these velocities are not modified 
when one makes the simplified hypothesis that the coefficient r\ is equal to zero. In [9] it was 
observed that the second sound velocity is much higher than that of the vortex waves, so that 
the small quantity r] should influence the two velocities (l5.14ll5Tl~5j) in a different way: negligible 
for the second sound velocity but relevant for the vortex waves. Regarding the attenuation 
coefficients (|5.16fJ5TTT|) . one sees that the first term in (I5.16P is identical to that obtained in 
|18j . when the vortices are considered fixed. The new term, proportional to A, comes from the 
interaction between second sound and vortex waves. 

It is to note that the first term (|5. 16[) produces an attenuation both to forward waves and 
to backward waves, while the second term contributes to the two kinds of waves in a opposite 
way, according to the sign of this term. Detailed measurements of the attenuation of second 
sound in directions parallel and orthogonal to the heat flux could allow us to establish the 
presence and the sign of this term. 

Note also that the second term of the dissipative coefficient kg is the same as the third 
term of fci 3 ' 4 \ but with an opposite sign. This means that this term contributes to the 
attenuation of the two waves in opposite ways; and its contribution depends also on whether 
the propagation of forward waves or of backward waves is considered. The first term of k s 3 ' 4 ^ 
produces always an attenuation of the wave, while the behavior of the third term is analogous 
to the first one. 



Second case: n orthogonal to q 

In order to make a more detailed comparison with the model studied in [U [9], one proceeds 
to analyze another situation, in which the direction of the wave propagation is perpendicular 
to the heat flux, that is, for example, assuming n = (0,0,1). This choice simplifies the 
system (15. 8p in the following form 

—ujpcyT + Kq 3 — ujpeiL = 
(— uj — iKBnyLa) qi — iKBnyaq\L = 
£KT - (w + \kB hv L(\ - a)) q 3 + X KL = 
iAL 3 / 2 ^ - (uj + irl 1 ) L + K J 3 = 
< rjKT + uKL + (-uj - v^kL) J 3 = (5.18) 

(— oj - ^kBhvL(1 - a)) q 2 = 
—i^\K,J\L + (—uj — zjikL) J\ = 
k (—uj - ijinL) J 2 = 

Note that, in contrast with what was seen before, but in agreement with the corresponding 
situation of the model described in [HE], here the transversal and the longitudinal modes in 
general do not evolve independently, as shown from the first five equations. However, one will 
see that this is the case if high-frequency waves are considered. 

As in the previous situation, assume that the values of the frequencies uj and of the real 
part of the wave number, k r , are high enough, such that the system (I5.18P may be simplified 
in the following form 
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—ujpcyf + k r q 3 — ujpeLL 
—uqi = 

£k r T — Loqz + x^rL = 

—ujL + k r J3 = 

r\k r T + i//c r L — uj J3 = 



(5.19) 







-cj<72 = — w J\ = — w J2 

Note that in this special case, as in the previous case and in [9], only the longitudinal modes 
are present, so that the dispersion relation assumes the form 



w w 



Vf + u 



V 



A 



12 



w 2 + V 2 2 v 



0, 



(5.20) 



pcv V An, 
which is similar to equation (|5.1ip . 

Now, the arguments suggested are the same to the previous situation, in fact, under the 
hypothesis (|5.12p . the dispersion relation (|5,20p takes the form 

w(w 2 - v){w 2 - Vi) = 0, (5.21) 

where V2 is the second sound velocity and yfv is the velocity of the vortex density waves 
in helium II. The conclusions which one achieves here are the same to those of the previous 
situation. Indeed, the modes corresponding to the solutions (|5.2ip are showed in Table [2j 
which, apart from the first column, are identical to those shown in Table [TJ 



w 


= 


to 1)2 


= ±v 2 


^3,4 


= ±v^ 




= V 


Qi = 





91 = 





f -- 

Q3 


= 
= 


f = 

93 = 


1> 

±V 2 pcvip 


f = 

93 = 


1 ( X-vP£l\ j, 
pc v { V'i-v ) V 

V 2 2 -u ^ 


L -- 


= 


L = 





L = 




J3 


= 


■h = 





Ji = 





Table 2: Modes corresponding to null velocity, second sound velocity, and vortex waves, respectively. 

Now, the same procedure than in the previous situation is followed, that is one supposes 
that all the quantities of the system (|5.18p . which don't appear in the system (|5.19p . are small 
enough compared to the other terms of the same system. Further, one also assume that the 
coefficient r] is not zero, but it has small enough values to be considered as a small perturbation 
to the physical system. Therefore, one assumes 



U) 



w 



— = w + 5 and K = k r + ik s , 



and substituting them in the dispersion relation of the system (15.18p . one finds the relation 
(I5.2ip . at the zeroth order in 5 and k s , and the following two expressions at the first order in 
5 and k q 



Wl,2 



n 



Ipcy ( W 



'1.2 



W 



3,1 



PtL 



X 



IV 



Wl,2, 



(5.22) 
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( 1 H r — ■ — j pe L - ] J tr :u . (o.23) 



and 




(5.24) 



(5.25) 



Q 2 s ) 

As regards the expression (|5.24p for the dissipative term kg , note that it is the same as 
the expression obtained when the vortices are assumed fixed [13, 20J, whereas the attenuation 
term /ci 3 ' 4 ^ is the same as the second term of /ci 3 ' 4 ^ of the first case (n parallel to qo). 



6 Conclusions 

The previous hydrodynamical model of inhomogeneous turbulent vortex tangle [H E] has been 
generalized, by adding the vortex flux J to the set of the independent variables E, q and L. 
In this new model, J is no longer described by a usual constitutive equation but it has its own 
dynamical equation. 

A set of evolution equations for E, q, L, and J subject to the restrictions of the second 
law of thermodynamics are studied and used to analyze the behavior of second sound and 
vortex waves, with special emphasis on their mutual couplings. This mathematical analysis 
may be useful in the interpretation of high-frequency second sound experiments, which play 
a key role as a tool for the probe of small spatial scales of turbulent vortex tangles. Indeed, 
the time derivative of J becomes relevant at high enough frequencies; in such a regime, the 
behavior of vortices becomes undulatory instead of being diffusive — the behavior assumed in 
the previous studies [H [9] . 

An interesting result is found from the comparison between the results of the wave propa- 
gation parallel and orthogonal to the heat flux. In fact, when the waves propagate orthogonal 
to the heat flux, the presence of the vortex tangle always causes an attenuation of the waves. 
But, when the propagation of the wave is collinear to the heat flux other terms are present. 
These terms have a positive or negative contribution depending on whether the direction of 
the wave is the same or opposite to the direction of the heat flux. 

Now, compare the results of the perturbed situation of the first case, when n is parallel to q, 
with those obtained in [HE]. The results of the comparison regarding the second case, n normal 
to q, are the same to those of the first case. As in [U [9], in this case one has the propagation 
of two kinds of waves, namely heat waves and vortex waves, which cannot be considered as 
propagating independently from each other. In fact, the uncoupled situation (equation (|5.13p ). 
in which the propagation of the second sound is not influenced by the fluctuations of the 
vortices, is no more the case when the quantities N\ = clkBhvL, N2 = nBuyaqi, N3 = AL 3 / 2 , 
N4 = JikJi, N5 = JikL, 1 and 77, appearing in the system (|5.9p . are considered. Indeed, 
from (|5.14fj5Tl~5"j) and from the results of [9] one makes in evidence that heat and vortex waves 
cannot be considered separately, that is as two different waves, but as two different features of 
the same phenomena. Of course, the results obtained here are more exhaustive than those of 
[HE]: in fact, comparing the velocities at the first order of approximation in both models, one 
deduces that the expressions (|5.14fj5TT5l) depend not only on the velocities of heat waves and 
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vortex waves, as in [U |9], but also on the coefficient rj, which comes from the equation (|5.1l i) 
of the vortex flux J, and whose physical meaning is a thermodiffusion coefficient of vortices. 
The fourth equation of the system (|5.9j) shows that the vortex flux J\ is not proportional to 
the heat flux, as it was assumed in [U [9] , but it satisfies an equation in which also the fields 
L and T, through 77, are present. 

It is to note that the attenuation of the second sound depends on the relative direction of 
the wave with respect to the heat flux: in some experiments this dependence was shown for 
parallel and orthogonal directions [21]. These results were explained assuming an anisotropy 
of the tangle of vortices, which corresponds to the presence of the parameter a in (|3.31j) . But, 
looking at the expressions (|5.16p and (|5.24p of the attenuation of the second sound in the 
high-frequency regime, one notes that these expressions are not equal after assuming a = 1/3 
(isotropy of the tangle). In particular, the term 



in ([5.16P causes a dependence of the attenuation depending on whether the wave direction 
agrees with the direction of the heat flux q or not. This term is absent if the wave propagates 
orthogonal to the heat flux. Note, in contrast, that the propagation speeds (|5.14j) and (|5.15j) 
for propagation direction n parallel to q coincide with (|5.22j) and (15.231) . respectively, for 
propagation direction n normal to q. Thus, the behavior of speed of propagation is isotropic 
and does not depend on the isotropy or anisotropy of the tangle. 

In conclusion, it could be that an anisotropy of the behavior of high-frequency second sound 
does not necessarily imply an actual anisotropy of the tangle in pure counterflow regime, but 
only a different behavior of the second sound due to the interaction with the vortex waves. 
This may be of interest if one wants to explore the degree of isotropy at small spatial scales. 
Of course, some more experiments are needed in order to establish the presence and the sign 
of these additional terms. 
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